Relaxation kinetics in two-dimensional structures 
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We have studied the approach to equilibrium of islands and pores in two dimensions. The two- 
regime scenario observed when islands evolve according to a set of particular rules, namely relaxation 
by steps at low temperature and smooth at high temperature, is generalized to a wide class of kinetic 
models and the two kinds of structures. Scaling laws for equilibration times are analytically derived 
and confirmed by kinetic Monte Carlo simulations. 

PACS numbers: 61.46.+W 
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I. INTRODUCTION 

A large amount of work has been done in the domain 
of atomic structure kinetics since Burton, Cabrera and 
Frank presented, in a seminal paper, the first serious at- 
tempt to model the behavior of atoms adsorbed onto a 
vicinal surface 1 . Nevertheless, several problems remain 
unsolved. Many important advances in the techniques 
for the observation and manipulation of atoms now allow 
dealing with structures of decreasing size and, with the 
help of powerful computation, gaining insight into the 
basic mechanisms governing atomic dynamics. (For a re- 
view, see for example Ref. yylj and references therein) . 

In the last few years, considerable efforts have been 
dedicated to understanding how nanostructures relax. A 
classical description of structure shape equilibration was 
developed by Herring, Nichols and Mullins (HNM) 5 . In 
this theory, the transport of mass responsible for relax- 
ation occurs in a smooth and continous way via migra- 
tion of border adatoms from regions of high curvature 
(high chemical potential fi) to regions of low curvature 
(low fi). It is assumed that structure border is rough 
enough to allow a continous (coarse-grained) description 
of it. However, as this is not valid for surfaces below the 
roughening temperature, the low temperature decay has 
been the subject of research both theoretically^*^ and 
exp er iment allyS*— 

Recently, studies of two-dimensional island shape re- 
laxation with a simple model has revealed new, unex- 
pected phenomena-^. It was shown that two qualita- 
tively different relaxation regimes exist. At high tem- 
perature, islands evolve toward equilibrium according to 
HNM. At low temperature, islands become faceted and 
a new driving mechanism appears. It was demonstrated 
that, in this case, shape relaxation occurs by steps, where 
the limiting process consists in the nucleation of new 
adatom rows on the flat edges of islands. In the above 
model adatoms lie on a triangular lattice and both equi- 
librium and kinetic properties depend on the single pa- 
rameter eo/ksT (where ks is the Boltzmann constant 
and T the temperature) since, for an adatom with rii 
nearest- neighbors (NN), rt^o is both the potential energy 
and the kinetic barrier controlling migration. 

This two-regime scenario is manifest in the depen- 



dence of the equilibration time t eq on temperature and 
island size N. At high temperature, t eq ~A 2 exp(3/3eo) 
(/3 = 1/fcuT), while at low T, t eq ~N exp(4/?e ). The 
regime to which relaxation of a given island belongs is 
determined by border roughness and depends not only 
on temperature but also on island size. In the T — N 
plane, the separation line between the two regimes is 
provided by the crossover island size N C (T) ~exp(/3e ), 
a rough indication of the size of the largest island that 
is fully faceted at T. Interestingly, if the temperature- 
dependent factors appearing in the leading terms of t eq 
are rewritten in terms of N c , the properly-scaled equili- 
bration time becomes a function of N/N c only, satisfying 
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A first problem of interest is the universality of this 
scaling law. The model described above is at best a first 
approximation to real situations and many details can 
be modified in order to improve it. For example, such 
a well-known effect as the Ehrlich-Schwoebel barrier- 
can be accounted for by introducing new kinetic barri- 
ers in a phenomenological way; in an even more accurate 
approach, minimum energy transition paths could be cal- 
culated from interatomic forces— As a different set of 
kinetic barriers will entail, to the least, new activation 
energies, the question remains whether the above scal- 
ing behavior is a consequence of the model or, to the 
contrary, would carry over to more general situations. 

Another closely-related problem is the relaxation of 
two-dimensional pores, i. e., islands of vacancies. In 
some sense, islands and pores are 'specular images' of 
one another. A given pore has the same boundary as the 
corresponding island, and so would be faceted or rough 
depending on whether the number of vacancies in it is 
smaller or larger than N C (T). It would be desirable to 
know if pore shape relaxation exhibits also a two-regime 
scenario and to analyze the possibility of scaling laws for 
the corresponding equilibration times. 

As far as we know, the questions above have not been 
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investigated, even in the case of simple models. In this 
work we study the main properties of shape relaxation 
in two-dimensional structures of adatoms. Our approach 
is two-fold: theoretical analysis and numerical simula- 
tions. In the analytic part, we treat the shape relaxation 
problem in a general way. Our study suggests the pres- 
ence of two regimes, regardless of the specific set of ki- 
netic barriers and, based on detailed-balance conditions, 
we show that similar scaling laws apply for both islands 
and pores. Since most of the arguments we invoke are 
heuristic, some independent confirmation of the derived 
properties is desirable. This is the goal of the numerical 
part. We analyze, using standard kinetic Monte Carlo 
(KMC) simulations, the relaxation of islands and pores 
according to two models considered in this paper, and 
compare the results with our theoretical predictions. 

The paper is organized as follows. In Sec.[n]we discuss 
the causes for scaling behavior and derive the asymptotic 
laws of scaling functions. The models we use are defined 
in Sec. Illll and the outcome of the KMC simulations are 
presented in Sec. IIVI Finally, we give our conclusions in 
Sec. E| 



to be created in order to attain the equilibrium shape; 
thus t eq ^N. 

The size exponents for equilibration time thus appear 
to be universal (1 in the FRM, 2 in the RRM). The ac- 
tivation energies however do depend in general on the 
kinetic barriers. The activation energy in the FRM cor- 
responds to the energy characteristic of row nucleation, 
while in the RRM it will be equivalent to the character- 
istic energy for diffusion along the island border. If we 
call these energies Ep and Ep for the FRM and RRM, 
respectively, the arguments presented so far allow us to 
write the ansatz 

r ^ e 0(E F+ e o ) for N_^ l 

te q ~\ 2 ■ (2) 

|^ ( e P(E R+ 2e ) for N_ ! 

At this point, it is clear that the occurrence of scaling 
is limited to the cases where 



II. ANALYTIC APPROACH 

In this section we investigate a generic model with 
adatoms lying on a triangular lattice where the total 
binding energy per NN pair is eo- Equilibrium properties 
depend on a single parameter, namely (3eo, but kinetics 
involve in general a greater number of them. Our aim 
here is to establish scaling laws for both island and pore 
relaxation times. 

Let us consider islands first. It is easy to see that many 
of the arguments presented in Ref. remain applicable 
regardless of the specific set of kinetic barriers. First of 
all, iV c (T), which separates rough and faceted islands at 
a given temperature, is the size at which the perimeter 
equals the average distance between border kinks. This 
does not depend on the transition barriers but only on 
the binding energies in the different configurations, so 
the form iV c ~ exp(/3eo) is still valid. Next, two qual- 
itatively different modes of relaxation should exist. On 
the one hand, for islands greater that N c , shape relax- 
ation dynamics is adequately described by the HNM the- 
ory. It consists of a set of equations for the temporal 
evolution of the border curvature, where the kinetics is 
dominated by perimeter diffusion and always leads to a 
t eq ^N 2 power law. We will call this regime rough re- 
laxation mode (RRM). On the other hand, mechanisms 
driving islands smaller than N c toward equilibrium arc 
also quite general. In this scenario, islands spend most 
of the time in highly-faceted configurations. We will call 
this regime faceted relaxation mode (FRM). Evolution oc- 
curs by steps consisting in nucleation and stabilization of 
new adatom rows. The time associated to each of these 
steps is proportional to facet length (^/V 1 / 2 ) and, given 
an initial configuration, a number ^N 1 / 2 of rows needs 



E F = E R + e ; (3) 

otherwise, the exponential factors in Eq. J5J) would be 
different and no natural time scale [such as ~7V 5 in Eq. 
||TJ], allowing scaled t eq to be expressed as a function of 
N/N c , would exist. Let us, then, analyze the connexion 
between the two energies. 
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FIG. 1: Schematic of the limiting processes for diffusion along 
a rough border. Adatom detaches from a kink and go into 
adjacent rows after jumping around a corner. All the sites 
below the full line, which represents the border, are occupied. 

As noted above, the limiting step in the RRM is 
adatom border diffusion. This is in turn limited by the el- 
ementary processes sketched in Fig. ^ hi which adatoms 
pass from one kink to another by jumping around a cor- 
ner; thus we can readily equal Ep and the energy of this 
process. 

For the FRM, the dynamics is determined by row nu- 
cleation, which occurs when two adatoms meet on a flat 
island edge. The nucleation rate may thus be calculated 
by the product of two factors: the rate of adatoms en- 
tering a facet and the probability that another adatom 
lies on the same facet. In this regime, the main sources 
of adatoms are facet ends, where the most weakly bound 
adatoms are found. Those adatoms move onto flat facets 
by basically the process sketched in Fig. ^ (A); the 
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first factor is thus ~ exp(— /3Er). The second factor is 
~ exp(— /3eo) because the system gains an energy eo when 
an adatom moves from a kink to a flat edge. Thus, we 
have Ef = Er + eo, i- e., Eq. J2J, and the scaling form 
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where the exponent a is 



Er 



(5) 



We consider now the case of pores. Since, as mentioned 
before, a pore will also be rough or faceted depending 
on whether it is larger or smaller than N c , it is worth 
analyzing each situation separately. 

For the relaxation of rough pores, the solution is im- 
mediate: the kinetic path is the same as that followed by 
the corresponding island (via a vacancy-atom exchange 
mechanism). This pore-island symmetry is clearly evi- 
dent in the HNM theory (see for example Ref. 0) , which 
involves only border curvature and perimeter atomic dif- 
fusion and hence does not depend on the kind of particles 
enclosed within the border. As a corollary, the equilibra- 
tion time for pores should also satisfy the asymptotic law 
established by the second equation in Eq. (J2J; this will 
be verified numerically in Sec. IIVI 



faceted islands, which again yields the power law t eq ~ N . 
In order to calculate the activation energy -which corre- 
sponds to the nucleation process illustrated in Fig. [2]- 
we consider first the inverse process, i. e., the migration 
of an adatom from site 2 to site 1. In this case, the 
adatom first leaves a kink, then diffuses along a flat bor- 
der and, finally, goes into an adjacent row by jumping 
around a corner. This is similar to the process sketched 
in Fig. ^ (B) and must therefore have an activation en- 
ergy Er. Knowledge of this last energy allows us to find 
that correponding to nucleation. They are related by the 
detailed-balance condition which, taking into accout that 
the binding energy difference between sites 1 and 2 is eo , 
finally leads to the result that the activation energies for 
pore shape relaxation are also connected by relation J3J. 

As an interesting consequence, the equilibration time 
for islands or pores, when properly scaled with iV", be- 
comes a function of N/N c only that satisfies the asymp- 
totic rules (QJ. We stress that even though the HNM 
theory implies an equivalent relaxation mode for both 
kinds of structures in the RRM, the same symmetry can- 
not be expected to hold in the FRM. The asymptotic 
forms give the exponents of the leading terms but dif- 
ferent prefactors may appear. In general, there will be 
one scaling function corresponding to islands and an- 
other scaling function for pores. For very small values 
of their arguments (N <C N c ) they will be parallel (on a 
log scale) but not necessarily equal, while at the other 
extreme (N^>N C ) the curves should collapse as required 
by HNM theory. Using KMC we will see, in Sec.[IV| that 
these properties of the scaling functions are verified. 




FIG. 2: Schematic of row nucleation in a pore. All adatoms 
are outside the line which represents the border of a fully 
faceted pore. 

For faceted pores, it is clear that relaxation will occur, 
as in the case of faceted islands, by steps consisting of 
the nucleation and stabilization of new rows. However, 
despite this qualitative analogy, the two processes are 
not equivalent but, rather, complementary; adatom mi- 
gration occurs from short to long facets in islands, while 
the opposite is true in pores. For the latter, row nucle- 
ation occurs when, as shown schematically in Fig. [21 an 
adatom that lies on a long facet moves to an internal 
corner of the pore. The new row will grow with the ar- 
rival of more adatoms from the now 'open' long facet, 
and become stable when fully filled. 

The similarity of the relaxation mechanisms allows us 
to apply, in the case of pores, the same kind of argu- 
ments we invoked for the dependence of t eq on the size of 



III. THE MODELS 

In order to asses the analytical predictions above, we 
investigated two different models, denoted I and II, which 
have the following common characteristics: The sub- 
strate is represented by a triangular lattice and the total 
binding energy per NN pair is eo. Adatom hoppings are 
only possible between NN sites and a given adatom in an 
initial site a jumps to a final empty site b with a prob- 
ability per unit time P a b = v exp(— f3E a i,), where v is a 
constant frequency and E a b the kinetic barrier. In or- 
der to avoid detachment, the jumps are forbidden if the 
number of NN adatoms in the final site is 0. As an addi- 
tional simplification, the motion of adatoms with initially 
more than 4 NN is also forbidden. This approximation is 
justified because of the high energy barriers of the corre- 
sponding processes; it is introduced in order to accelerate 
the simulations. Differences between models I and II lie 
only in the set of kinetic barriers. 

For Model I, which is the same as the one used in Ref. 
00, E a b depends only on the number n a of NN adatoms 
in the initial site a. It takes the value E a t, = n a £o when 
1 < n a < 5 and E a b = e/10 when n a = l, regardless of the 
number of NN adatoms in b. 

For Model II, we use a more accurate set of activation 
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TABLE I: Kinetic barrier for Model I (AE 7 ) and Model II 
(AE 11 ) as a function of the number of NN before (m) and 
after (n/) the jump. For Model I, the barriers depend only 
on rii. Detachments, as well as jumps of adatoms with rii >4, 
are forbidden in both models. 



IV. NUMERICAL RESULTS 



FIG. 3: The energy landscape for Model II. 



energies, obtained from a detailed study of the kinetic 
barriers for adatoms lying on a triangular lattice and in- 
teracting via a Lennard- Jones potential^. It is useful to 
define this model with the help of Fig- El where the differ- 
ent barriers are represented. In the top part of the figure, 
adatoms are represented by circles while crosses indicate 
possible binding sites. The bottom part shows the corre- 
sponding energy landscape, i. e., kinetic barriers separat- 
ing adjacent local minima. Thus, in this model, we take 
into account the final-state configuration when evaluat- 
ing the transition-state energy. (Note however that the 
height of the barrier depends on the number of NN in the 
initial and final states but not on their precise position) . 
Model II thus incorporates such important features as 
the Ehrlich-Schwoebel barrier. This is apparent if we 
look at Fig. |3J the barrier from 2 to 1 is lower than the 
barrier from 2 to 3; the barrier from 4 to 3 is higher 
than the barrier from 4 to 5, etc. The final parameters 
in this model are i?i = 0, £2=0. 5eo, E^^eo, £4 = 1. 5eo, 
E$ = 2eo, i?6=4eo- In order to facilitate comparison, Ta- 
ble lists, for both models, the kinetic barriers corre- 
sponding to a given elementary jump as a function of the 
number of NN before (n^) and after (nj) the hopping. 

Before closing this section, we calculate the activation 
energies for each model as derived from the analysis of 
Sec. m Remember that En is the energy characteristic 
of the process sketched in Fig. ^ (A). For Model I, this 
process corresponds to an adatom with 3 NN before the 
jump, so E R — 3eo, and, using Eqs. J5J and ©, Ep=4:Co 
and a 1 = 5, as obtained in Ref. |6j|7J (but not yet tested 
with pores). For Model II, the process in question is (for 
example) the transition from 4 to 2 in Fig. [3 which leads 
to EK = 2e , £j7 = 3e and a 11 =4. 



As mentioned in Sec.[fl a numerical confirmation of the 
results of Sec.[nl obtained in a heuristic manner, is in or- 
der. In this section, we do this using KMC simulations 
to explore shape relaxation according to models I and 
II. We consider islands and pores with sizes between 400 
and 6500 particles and perform, for each model, standard 
KMC simulations. The values of j3 (in 1/eo units) range 
between 2 and 17. The results shown in this section cor- 
respond to averages over a number of samples, varying 
from 4 for the largest structures to 30 for the smallest 
ones. 

Following Ref. 00, we use the aspect ratio a, defined as 
the ratio of the x and y gyration radii, to characterize the 
state of the structures. We start each simulation with an 
aspect ratio around 10 and define the equilibration time 
t eq as the time at which a first becomes less than 1. First 
of all, we analyze the dependence of t eq on N (for fixed j3) 
and, separately, on (3 (for fixed N). Later in this section, 
we will check the scaling laws. 

Let us start with Model I. In order to test the behav- 
ior of teq as a function of TV, we plot in Fig. 0] (o) \n.t eq 
against In N for several temperatures. Filled and empty 
symbols correspond to islands and pores, respectively. At 
high (low) enough temperature, t eq is expected to scale 
as A^ 2 (N). This scaling behavior is represented in Fig. 

01 (a) by the lines with slopes of 2 (lower) and 1 (up- 
per), which bracket the KMC results corresponding to 
the highest and lowest temperatures used in our simu- 
lations, respectively. Using these lines as a guide, it is 
indeed clear that the the size exponent changes from 1 to 

2 when the temperature is increased, consistent with the 
numerical findings. Note that the agreement is as good 
for pores as it is for islands. 

We check now our predictions for the activation ener- 
gies. To do so, in Fig. 01 (6), h\t eq is plotted against /3 
for different sizes. The asymptotic behaviors are repre- 
sented by the straight lines which have slopes of 3eo (low 
(3) and 4e (high (3). The KMC data again show excel- 
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FIG. 4: Equilibration time for Model I. (a) t eq as a function of 
JV for f3 = 3.9 (circles), 5.8 (pentagons), 11.6 (down-triangles), 
12.9 (squares), 14.5 (diamonds), 16.6 (up-triangles); the lines 
have slopes of 1 (upper) and 2 (lower), (b) t eq against f3 for 
N = 490 (squares), 1000 (triangles), 4000 (circles); the lines 
have slopes of 3 (left) and 4 (right). Filled symbols correspond 
to islands and empty ones to pores. 



FIG. 5: Equilibration time for Model II. (a) t eq as a 
function of N for /3 = 3.9 (circles), 4.6 (pentagons), 5.8 
(down-triangles), 10.6 (squares), 11.6 (diamonds), 12.2 (up- 
triangles); the lines have slopes of 1 (upper) and 2 (lower), 
(b) teq against /3 for N — 490 (squares), 1000 (triangles), 
4000 (circles); the lines have slopes of 2 (left) and 3 (right). 
Filled symbols correspond to islands and empty ones to pores. 



lent agreement with the predicted behavior, to within a 
trivial prefactor (i. e., slopes agree). Let us remark that 
the results for island relaxation are given here to allow 
comparison with the (new) results for pores. We note 
however a small discrepancy between our numerical re- 
sults and those of Ref. 0: although the size exponents 
and activation energies coincide, our equilibration times 
are found to be lower (by a factor of ~ 10) than those 
reported in Ref. 

The corresponding results for Model II are presented in 
Fig. El Note that, also in this case, the theory-simulation 
agreement is good for both types of structures. Thus, 
Figs. 01 and [S] confirm our predictions that, on the one 
hand, the size exponents are universal (1 and 2) and, on 
the other hand, En and Ep depend on the specific set of 
microscopic kinetic barriers. 

Finally, in order to test the specific scaling laws, we 
plot, on a log-log scale, (t eq /Nc 5 ) against (N/Nc) for 
Model I (Fig.EJl and {t eq /Nc 4 ) against {N/Nc) for Model 
II (Fig.0). We have used a 1 and a 11 as calculated at the 
end of Sec. IIIII According to the analysis presented in 
Sec. 1^ teq/Nc should scale with exponents 1 and 2 for 
very small and very large values of its argument (N/N c ). 
In Figs. El and [7| we have drawn lines with slopes of 1 and 
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FIG. 6: Model I scaling of log 10 (t eq /N^) against log 10 (JV/JV c ), 
for islands (filled symbols) and pores (empty symbols). The 
lines have slopes of 1 and 2. 



2 to indicate these asymptotic behaviors, which are evi- 
dently closely verified. For concreteness, in these plots, 
we use the same N c = .25 exp(/3eo) as in Ref. |(| the exact 
value of the geometrical prefactor should not affect the 
final result in a significant manner. 

It is interesting to note that, besides confirming the 
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FIG. 7: Model II collapse of data for both islands (filled sym- 
bols) and pores (empty symbols). The lines have slopes of 1 
and 2 



expected asymptotic power laws, in both figures all the 
points corresponding to islands collapse on one curve, 
and all the points corresponding to pores collapse on an- 
other one (depending via a on the microscopic details of 
the model), thus giving additional support to the scaling 
ansatz. Furthermore, although for each model the island 
and pore scaling functions are different in the FRM, they 
become equivalent in the RRM, as required by the perfect 
pore-island symmetry of relaxation in the latter regime. 



V. CONCLUSIONS 

By considering a generic model of adatoms lying on 
a triangular lattice, we have studied the problem of 
shape relaxation of islands and pores. The arguments 
employed in this work allows the main properties of 
the equilibration time t eq to be calculated as a func- 
tion of temperature and size. Because our arguments 



are somewhat heuristic, KMC simulations, using two 
different kinetic models, were also carried out. The 
numerical results confirm our theoretical predictions. 
For both islands and pores two qualitatively different 
modes of relaxation (FRM and RRM) are found, as 
well as the line N C (T) that separates these regimes 
in the temperature-size plane. It was shown that, 
although size exponents are universal (1 in FRM, 2 
in RRM), the activation energies corresponding to 
each mode depend on the microscopic details of the 
kinetic model. Scaling behaviors were found for the 
equilibration time: the properly scaled t eq depends also, 
via a, on the set of elementary kinetic barriers of the 
model. Furthermore, for a given model, the specific 
scaling function for islands is in general different from 
the one corresponding to pores. While in the FRM 
both functions scale with the same exponent (= 1) as 
a consequence of the detailed-balance condition, in the 
RRM the two functions collapse because of the perfect 
island-pore symmetry of HNM theory^. 
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